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Nonhomogeneous boundary conditions for the 

SPECTRAL FRACTIONAL LAPLACIAN 

NICOLA ABATANGELO, LOUIS DUPAIGNE 


Abstract 

We present a construction of harmonic functions on bounded domains for the spectral fractional 
Laplacian operator and we classify them in terms of their divergent profile at the boundary. This 
is used to establish and solve boundary value problems associated with nonhomogeneous boundary 
conditions. We provide a weak-L^ theory to show how problems with measure data at the boundary 
and inside the domain are well-posed. We study linear and semilinear problems, performing a sub- 
and supersolution method, and we finally show the existence of large solutions for some power-like 
nonlinearities. 


1 Introduction 

Given a bounded domain LI of the Euclidean space , the spectral fractional Laplacian operator 
s G (0,1), is classically defined as a fractional power of the Laplacian with homogeneous 
Dirichlet boundary conditions, seen as a self-adjoint operator in the Lebesgue space L^(n), see (2) be¬ 
low. This provides a nonlocal operator of elliptic type with homogeneous boundary conditions. Recent 
bibliography on this operator can be found e.g. in Bonforte, Sire and Vazquez [5], Grubb [15], Gaffarelli 
and Stinga [7], Servadei and Valdinoci [18]. 

One aspect of the theory is however left unanswered: the formulation of natural nonhomogeneous 
boundary conditions. A first attempt can be found in the work of Dhifli, Maagli and Zribi [11] . The 
investigations that have resulted in the present paper turn out, we hope, to shed some further light on 
this question. We provide a weak formulation, which is well-posed in the sense of Hadamard, for linear 
problems of the form 

[ (-^ln)"w= M in fl, 

I ^ = C on afl, 

where hi is a reference function, see (8) below, with prescribed singular behaviour at the boundary. 
Namely, hi is bounded above and below by constant multiples of where 

S{x) := dist(a:, dLl) 

is the distance to the boundary and the left-hand side of the boundary condition must be understood as a 
limit as S converges to zero. In other words, unlike the classical Dirichlet problem for the Laplace operator, 
nonhomogeneous boundary conditions must be singular. In addition, if the data p,, are smooth, the 
solution blows up at the fixed rate . This is very similar indeed to the theory of nonhomogeneous 

boundary conditions for the classical (sometimes called "restricted") fractional Laplacian - although in 
that case the blow-up rate is of order - as analysed from different perspectives by Grubb [14] and 

the first author [1] of the present note. In fact, for the special case of positive s-harmonic functions, that 
is when ^ = 0, the singular boundary condition was already identified in previous works emphasizing the 
probabilistic and potential theoretic aspects of the problem: see e.g. Bogdan, Byczkowski, Kulczycki, 
Ryznar, Song and Vondracek [4] for an explicit example in the framework of the classical fractional 
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Laplacian, as well as Song and Vondracek [21], Glover, Pop-Stojanovic, Rao, Sikic, Song and Vondracek 
[12] and Song [20] for the spectral fractional Laplacian. 

Turning to nonlinear problems, even more singular boundary conditions arise: in the above system, 
ii fi = —u^ for suitable values of p, one may choose C = +oo, in the sense that the solution u will blow 
up at a higher rate with respect to controlled by the (scale-invariant) one Note 

that the value C = -l-oo is not admissible for linear problems. This was already observed by the first 
author in the context of the fractional Laplacian, see [2], and this is what we prove here for the spectral 
fractional Laplacian. Interestingly, the range of admissible exponents p is different according to which 
operator one works with. 

Main results 

For clarity, we list here the definitions and statements that we use, with reference to the sections of the 
paper where the proofs can be found. First recall the definition of the spectral fractional Laplacian: 

Definition 1. Let C a bounded domain and let be a Hilbert basis of L^{H) consisting 

of eigenfunctions of the Dirichlet Laplacian —Aj^, associated to the eigenvalues Xj, j € N, i.ef (pj G 
iLp (n) n C°°{Lt) and —L\ipj = Xjipj in Q. Given s G (0,1), consider the Hilbert space^ 

{ OO OO 

i=i j=o 

The spectral fractional Laplacian of u G H(2s) is the function belonging to Lf{Ll) given by the formula 

OO 

(-A|o)®w = (2) 

i=i 

Note that C H(2s) ^ Lf{Vl). So, the operator (—Ajj^)* is unbounded, densely defined and 

with bounded inverse (—A[q)~^ in L‘^{Lt). Alternatively, for almost every x G LI, 

= Py [ [u{x) - uiy)]J{x,y) dy + k{x)u{x), (3) 

Jn 

where, letting p(i(t,x,y) denote the heat kernel of — Ajj^, 

= = f(i^<“> 

are respectively the jumping kernel and the killing measure,^ see Song, Vondracek [21, formulas (3.3) 
and (3.4)[. For the reader’s convenience, we provide a proof of (3) in the Appendix. We assume from 
now on that 

n is of class . 

In particular, sharp bounds are known for the heat kernel po, see (27) below, and provide in turn sharp 
estimates for J{x,y), see (33) below, so that the right-hand side of (3) remains well-defined for every 
X gLI under the assumption that u G fl L^{fl, 6{x)dx) for some e > 0. This allows us to define 

the spectral fractional Laplacian of functions which do not vanish on the boundary of II. As a simple 
example, observe that the function u = 1 does not belong to H(2s) if s > 1/4, yet it solves (1) for p, = «: 
and C = 0- 

^See Brezis [6, Theorem 9.31]. 

■^When is smooth, H[2s) coincides with the Sobolev space if 0 < s < 1/4, Hqq{Q.) if s E {1/4, 3/4}, 

otherwise; see Lions and Magenes |16, Theorems 11.6 and 11.7 pp. 70-72]. 

®in the language of potential theory of killed stochastic processes. Note that the integral in (3) must be understood in 
the sense of principal values. To see this, look at (33). 
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Definition 2. The Green function and the Poisson kernel of the spectral fractional Laplacian are defined 
respectively by 

Ghix,y) = pn{t,x,y)G~^ dt, x, y G n,x y, s € {0,1], (5) 

r(s) Jo 

where pn denotes the heat kernel of — A|q, and by 

PQix,y) ■■= -Gf^{x,y), xGn,yGdn. (6) 

where v is the outward unit normal to dVl. 


In Section 2, we shall prove that Pq is well-defined (see Lemma 12) and review some nseful identities 
involving the Green function Gq and the Poisson kernel Pq. Now, let ns define weak solntions of (1). 

Definition 3. Gonsider the test function space 



T{n) := (-A|^)- 

~^CT{Pl) 

(7) 

and the weight 





hi{x)= PQ{x,y)da{y), x G Pt. 

(8) 


JdQ 



Given two Radon measures 

y G A4{Pl) and ( G JA{dPl) with 



/ S{x) dlyl{x) < oo, 

Jq 

ICl(afi) < oo. 

(9) 

a function u G Lj^^{Pl) is a 

weak solution to 




r(-A|^)^u= y 

in PI 



\ 

on dPl 

(10) 


if for any 'if gT{P i), 

[ w(-A|o)> = / 'ifdfi- f ^ df. (11) 

Jq Jn Jon 

We shall prove that T{Pi) C (^^(n), see Lemma 26, so that all integrals above are well-defined. 
Eqnation (11) is indeed a weak formnlation of (10), as the following lemma shows. 

Lemma 4. 1. (weak solutions are distributional solutions) Assume that u G L(^^{fl) is a weak solution 

of (10). Then in fact, u G L^{Pl, 6{x)dx) and (—Aj^)^ 'u = pL in the sense of distributions i.e. for 

any %f G Cf°{Pl), -—is bounded and 

[ u{-A\Q)"if= f ifdy, 

Jfi Jfi 


moreover the boundary condition holds in the sense that for every (f G C{Pl) 


lim — 
40 t 


'{S(x)<t} 


u{x) 

hi{x) 


(f{x) da{x) 



<f{x) dC{x) 


whenever y, G A4(n) satisfies (9) and ( G L^{dPl). 

2. (for smooth data, weak solutions are classical) Assume that u G L\^JPi) is a weak solution of (10), 
where y G C°‘{Pl) for some a such that a-|-2s ^ N and ( G C{dPl). Then, (—A|q)*u is well-defined 
by (3) for every x G PI, (—A|Q)*M(a;) = y{x) for all x G PI and for all xq G dPl, 


X 


u{x) 

lim 

^ xo hi{x) 


({xo)- 
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3. (classical solutions are weak solutions) Assume that u G that u/hi G C(0). Let 

fx = (—A|q)*u be given by (3) and ( = Then, u is a weak solution of (10). 

We present some facts about harmonic functions in Section 3 with an eye kept on their singular 
boundary trace in Section 4. 

We prove the well-posedness of (10) in Section 5, namely 

Theorem 5. Given two Radon measures /i G A4(n) and ( G Ai{dVL) such that (9) holds, there exists a 
unique function u G which is a weak solution to (10). Moreover, for a.e. x gLI, 


l{x) = [ Gf^{x,y) dfj.{y) + [ PQ{x,y) dC{y) 
Jq JdQ. 


and 


II ^ II (f2,(5(a;)(ia:) — s) (II^MllM(n) + IICIlA<(an)) • 

In addition, the following estimates hold. 


|w||LP(n,5(x)dx) < C'i||^/i||x(n) 


if ( = 0 and p G 


1 , 


iV + 1 


( 12 ) 

(13) 

(14) 

(15) 

(16) 


_ ’ iV + 1 - 2s 

||m||c“(u) < (^2 (||m||l=<=(c^) + ||CllAi( 9 n)) ifujCCQ.andaG{ 0 , 2 s) 

I|u||c2»+“(:j) < C'a (||mI|c“(u) + IICIl7W(9a)) if uj CC It and 2s + a 

In the above Ci = Ci (O, N, s,p), C2 = C'2(n, w, N, s, a), C3 = C3{Il, ut, N, s, a). 

In Section 6 we solve nonlinear Dirichlet problems, by proving 

Theorem 6. Let g(x, t) : H x R’*' —> R+ be a Garatheodory function and h : R+ —^ R+ a nondecreasing 
function such that g[x,0) = 0 and for a.e. x G It and all t > 0, 

H If g(x,t) < h{t) where h{S~^‘^~'^^'^)6 G L^ (It) 

(-^ln)"«= -9{x,u) 


Then, problem 


hi 


= C 


in n 
on dit 


(17) 


has a solution u G L^{It, 5(x)dx) for any f G C(dlt),( > 0. In addition, if t g(x,t) is nondecreasing 
then the solution is unique. 


Finally, with Section 7 we prove 
Theorem 7. Let 


pG 1 + s, 


1 — s 


Then, there exists a function u G L^(It,5{x)dx) n C°°{It) solving 

(—A|Q)®rt= —u^ in It, 

— = +00 on dIt 

hi 


(18) 


in the following sense: the first equality holds pointwise and in the sense of distributions, the boundary 
condition is understood as a pointwise limit. In addition, there exists a constant C = C{It, N, s,p) such 
that 

0 < u < CS~ p-i. 

Remark 8. If one replaces the spectral fractional Laplacian by the classical fractional Laplacian (—A)®, 
then an analogous of Theorem 7 holds for 

1 


p G 1 + 2s, 


1 — s 


see [2, Theorem 3 and Example 22]. We suspect that both exponent ranges are optimal. 
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2 Green function and Poisson kernel 

In the following three lemmas"'^, we establish some nsefnl identities for the Green function defined by (5). 
Lemma 9. ([13, formula (17)]) Let f G L‘^{Lt). For almost every x € fl, Gfi{x, •)/ G L^iVi) and 

(-^lo)”V(a:) = / G^^{x,y)f{y) dy fora.e.xeVt. 

Jn 

Proof. If ipj is an eigenfunction of — A|j^, then 

/ Gf^{x,y)ipj{y) dy = / / pn{t, x,y)ipj{y) dy dt 

Jn Jo ^ [s) Jn 


e ^^*ipj{x) dt = ^ 


roo ^s-1 

lo W) 


r(s) 


Pjix) F e dt = A-Vj(a;) = (-Ain) Pjix) 


By linearity, if / G L^(0) is a linear combination of eigenvectors / = J2j=i fjPP then 


and so, letting 


we have 


/ Gl^{x,y)^fjipj{y) dy = ^ijA-Vj(a;). 

<^hf ■■= [ GU;y)f{y)dy, 

Jn 


M 


i=i 


lL2(n)- 


(19) 


( 20 ) 


Thus the map Gq : / i—^ Gr^f uniquely extends to a linear isometry from L^{Ll) to H{2s), which 
coincides with (—A|n)~*. It remains to prove that the identity (19) remains valid a.e. for / G Lf{Ll). 
By standard parabolic theory, the function (t,x) i—>■ J^pQ{t.,x,y)dy is bounded (by 1) and smooth in 
[0,r] X u) for every T > 0, oj CC LI. Hence, for every x G LI, Gn(x, •) G L^{Ll). Assume first that 
f = tp G G)?°(H) and take a sequence in the linear span of the eigenvectors such that 

{V’fcIfeeN converges to ip in L^iLt). The convergence is in fact uniform and so (19) holds for f = ip- Indeed, 
by standard elliptic regularity, there exist constants G = G{N, Ll),k = k{N) such that any eigenfunction 
satisfies 

W^PjWL-m < (CA,)'=||(^,|U.(n) = (CA,)^ 

In particular, taking G larger if needed. 


S 


L“(n) 


< {gXjV- 


( 21 ) 


Now write ip = '^jPj ^cid fix m G N. Integrating by parts m times yields 


iPj = / iPpj = - 

Jn 


Xj Jn 


■0Av^j = / Aipipj = ...= 


{-ly 


•j Jn 


A™'!/' Fj 


which implies that 


I 0 il < 


l|A'"0|U.(o) 


\ m 

A- 


( 22 ) 


i.e. the spectral coefficients of ip converge to 0 faster than any polynomial. This and (21) imply that 
{0fc}fcgN converges to ip uniformly, as claimed. 


Take at last f G and a sequence {fk} 


fceN 


G)?“(0) of nonnegative functions such that {fk} 


fceN 


converges to |/| a.e. and in L^(H). By (20), ||GQ/fc||i 2 < ||/fc||i 2 and by Patou’s lemma, we deduce that 
Gq{x, •)/ G L^{Lt) for a.e. x G LI and the desired identity follows. □ 


which hold even if the domain £7 is not C^J. 
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Lemma 10. ([13, formula (8)]) For a.e. x,y € 

[ Gl,-%x,OGhi^,y)d^ = Gh{x,y) (23) 

JQ, 

Proof. Clearly, given an eigenfunction ipj, 

(-AIJ-^(-AI^)^-V, = A-^Af V, = (-A|n)-V, 

so (—A|j^)~® o (—= (—A|q)~^ in L^(VL). By the previous lemma and Fubini’s theorem, we 
deduce that for tp £ Lf(yi) and a.e. x £Vl, 


[ "{x,0Gn{^,y)p{y) dfdy= f G\^{x,y)p{y) dy 

Jn^ Jn 

and so (23) holds almost everywhere. □ 

Lemma 11. For any ip € G^{n), 

(-A|J> = (-A)o(-AIJ^-V = (-A|J*-^o(-A)^ (24) 


Proof. The identity clearly holds if ip is an eigenfunction, li ip £ (7“ (11), its spectral coefficients have 
fast decay and the result follows by writing the spectral decomposition of ip. Indeed, thanks to (22) and 
(21), we may easily work by density to establish (24). □ 

Let us turn to the definition and properties of the Poisson kernel. Recall that, for x £ fl,y £ dfl, the 
Poisson kernel of the Dirichlet Laplacian is given by 

d 


Lemma 12. The function 


Phi.x,y) = - — G}^{x,y). 
dvy 


Pa{x,y) ■■= —^Glt{x,y) 


is well-defined for x £ fl,y £ dfl and Pq(x.,-) £ G{dn) for any x £ fl. Furthermore, there exists a 
constant G > 0 depending on N, s, fl only such that 


1 


d{x) 


G\x-y\N+2-2s 


< PQix,y) < G- 


S{x) 


_ ryi\N-\-2 -2S 


x-y\ 


and 


[ Gl,-%x,0Pni^,y)df = Phix,y). 
Jn 


(25) 


(26) 


Remark 13. When H is merely Lipschitz, one must work with the Martin kernel in place of the Poisson 
kernel, see [12]. 

Proof, of Lemma 12. Take x, z £ Ll, y £ dLl. Then, 

pn{t,x,z)t ^dt= pn{\z-x\^T,x,z)T dr 

S{z) r{s)S{z)Jo T{s)S{z)Jo 

Since fl has boundary, given x £ ft, pq{-,x,-) £ C^((0,+oo) x 11) and the following heat kernel 
bound holds (cf. Davies, Simon and Zhang [9,10,24]): 


Hx)6{y) 

t 


A 1 


__^e 1"^ <pn(t,x,y) < 


d{x)6{y) 


A 1 


_£2_p-C2\x-y\^/t 

tNl2 


e (27) 


where ci,C 2 are constants depending on 11, iV only. So, 


(5(z) 


-pn{\z - x\'^t,x,z)t’^ ^<C|2;- 


s-l ^ r<\^ ™|2s-Ar-2jj-^^^s-2-Af/2g-02/T 


(28) 
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and the reverse inequality 


^ I |2s—A^ —2^/ N S—2 — N/2 —IKcot) ^ ^1 /i |2 n s—1 

-^\z-x\ 5{x)t ' e — ——p^[\z-x\t,x,z)t 

also holds for t > S{x)5{z)\z — x\~^. As z ^ y € dil, the right-hand-side of (28) obviously converges in 
L^(0, -fcx), dr) so we may apply the generalized dominated convergence theorem to deduce that Pq{x, y) 
is well-defined, satisfies (25) and 


Pa{x,y) 


-^Gh(.x,y) = lim 


diy. 


z^y 


Gq{x,z) 

6{z) 


1 

W) 


f°° d 

/ -7^Pni.t,x,y)P~^ dt. 

Jo 


From this last formula we deduce also that, for any fixed x &Q,, the function Pq{x, •) € C{dil): indeed, 
having chosen a sequence {yk}k^fi C dfl converging to some y G dfl, we have 


\PQix,yk) - Pnix,y)\ < 


m 


d . . d . , 

T^pn[t, X, yk) - -^pn{t, x, y) 


dv, 


duy 


P ^ dt 


where, by(27) 


dv. 


-Pn{t,x,y) 


Cl S{x) C 2 \x-y\'^/t ^ Cl 6{x) ^.^sjxPlt 

- ■tN/2+1 ^ - tN/2+1 ^ 


for any y G dfl, 


so that |Fn(a;, yk) — PQ(a;, ?/)| —>• 0 as fc t oo by dominated convergence. 

_2 

By similar arguments, is a continuous function on 12 \ {(a:, y) : x = y). And so, by (23), we have 

- 7 ^ / G\^\x,OGh{U)d^= P^{x,y) 

OVy Jq 

Let us compute the derivative of the left-hand side alternatively. We have 






(C,2) 


S{z) 


d^ = 


/R+ xH 


f{t,^,z) dtd^, 


where, having fixed a: G 12, 


pn(t,^,z)P <G|x-C| 


2s-N 


d{x)S{0 




A 1 


r(s)(5(z) 

For fixed e > 0, and ^ G 12 \ B{y, e), z G B{y, e/2), we deduce that 

f{t,^,z) < C\x-^f^-^p-^-^Ge-<^e/t g ii((o,+oo) X 12) 

Similarly, if 2 > e. 




f{t,^,z)<C\x-^\'^^ ^ G L^((e,-hcx3) X 12) 


Now, 


p-^-N-2e-‘^-^^^dt < 1 ^ - 


^+oo 


Hence, there exists a constant G > 0 independent of e such that 

[ f{t,^,z)dtd^<Ce^+^^ 

J {0,e)xB{y,e) 

It follows from the above estimates and dominated convergence that 


Ph{x,y) = lirn ( f{t,^,z)dtd^= [ G}^ ^{x,^)P^{^,y)d^ 

J (o,+oo)xn Jn 


i.e. (26) holds. 


( 29 ) 

□ 
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Remark 14. Thanks to the heat kernel bound (27), the following estimate also holds. 


C \x-y\ 


N-2s 


1 A 




-yp / 


< 


c 


\x-y\ 


N-2s 


6{x)6{y) \ 

\x-y\^)- 


( 30 ) 


for some constant C = C{fl, N, s). Also observe for computational convenience that 


1 


1 A 


Hx)S{y) \ ^ S{x)S{y) < Pi a 

k - yP / S(x) ^(y) + k - yp “ V 1^; - yP J ' 


3 Harmonic functions and interior regularity 

Definition 15. A function h G L^{fl,S{x)dx) is s-harmonic in fl if for any if G there holds 

[ ^(-A|n)> = 0. 

Jn 

The above definition makes sense thanks to the following lemma. 

Lemma 16. For any ip G C'“(fl), (—A|q)^'!/; G C'o(n) and there exists a constant C = C{s, N,il,ifj) > 0 
such that 

|(-A|j.j)>| < C5 inn. (31) 

In addition, if tp > 0, tp ^ 0, then 


(—AIq)*!/) < —CS inn\supp'tp. 
Proof. Thanks to (21) and (22), (—A|q)*'!/; G C'o(n) and 

(-A|o)> 


(32) 


i=i 


fl 

6 


L°“(n) 


< oo 


and (31) follows. Let us turn to the case where ip > 0, pj ^ 0. By the heat kernel bound (27), there 
exists C = C{fl,N,s) > 0 such that 


C\x-y\ 


N+2s 


iFlihpAi 

\x-y\ 


< J{x,y) < 


C 


\x-y\ 


N+2s 


S{x)S{y) 

|a; -yP 


A 1 


(33) 


Now, we apply formula (3) and assume that x G n\ supp-0. Denote by x* a point of maximum of ip and 
let 2r = dist(a;*,supp'0). Then for y G Br{x*), it holds ip{y), S{y) > ci > 0, |a; — ?/| < C 2 and so 


(-A|n)Xa;) = - [ i’iy)J{x,y) dy 

Jn 

i’i.y) 


< -c 

JBr{x>) \x - y\ 

< -CS{x). 


iV+ 2 s 


Hx)Siy) 
|a; -yP 


A 1 


dy 


□ 


Lemma 17. The function Pq{-,z) G L^{n,S{x)dx) is s-harmonic in D for any fixed z G dil. 
Proof. Thanks to (25), Pq{-,z) G L^(fl, P(x)dx). Pick ip G Cp!°{n) and exploit (24): 


[ P^(.,z)(-A|J>= / P^(.,z)(-Ay*-^[-AV^]. 
Jn Jn 


/Q, J Q, 

Applying Lemma 9, the Fubini’s Theorem and (26), the above quantity is equal to 


/ -Fn(-,^:) (-A)V'= 0. 
Jn 


□ 




























Lemma 18. For any finite Radon measure G let 


h{x) = [ P^{x,z)dCiz), xGft. 

Jan 


(34) 


Then, h is s-harmonic in 12. In addition, there exists a constant C = C{N, s,Tl) > 0 such that 

\\h\\L'^(Q^5{x)dx) < C'IICIIx(aa)- (35) 

Conversely, for any s-harmonic function h € 5(x)dx), h > 0, there exists a finite Radon measure 

C G M{dQ), C > 0> such that (34) holds. 

Proof. Since P{x, •) is continuous, h is well-defined. By (25), 

d\C\{z) 


S{x)\h{x)\ < C 


N-2s 


Jan \x-z, 

so that h € L^{Tl,5{x)dx) and (35) holds. Pick now fj € 

[ hix)i-A\^y'ilj{x) dx = [ ([ P^{x,z){-A\^yilj{x) dx\ dC{z) = 0 
Jn Jan \Jn ) 

in view of Lemma 17. Conversely, let h denote a nonnegative s-harmonic function. By Definition 15 and 
by equation (24), we have for any ip € C“(D) 

0= / h{x)i-A\^YiP{x) dx = f h{x) {-A\f^Y~^ o {-A)iPix) dx = 

Jn Jn 

Gh~"ix,^)h{0df] {-Aip{x)dx, 


n \Jn 


so that JqGq ^(x,^)h(^)df is a (standard) nonnegative harmonic function. In particular (cf. e.g. [3, 
Corollary 6.15]), there exists a finite Radon measure ( € A4(dfl) such that J^Glfi^{x,£,)h{f^)df = 
fan y)^C(y)- We now exploit equation (26) to deduce that 


f gIyYaoIhO- f P^(Cy)d((y) 

Jn L Jan 


df = Q. 


Since 


m 

ll —S 


ifi(x) [f Gli dx= ( ^(-Als.;)'’ Vi = Ti^ / hifi 

\Jn J Jn Jii 


< oo. 


it holds JqGq ^(x,^)h(^)d^ G C^(fl) fl L^(fl,S(x)dx). Thanks to (30), we are allowed to let Gq act on 
it. By (23), this leads to 


Gh(A 0 [/^(e) - [ y) rfc(2/)l = 0. 

! L J 


df 


Take at last ip G Gf°(I2) and = (—Ajj.^) ^ip. Then, 

0= / (fix) [ Gl^{x,^) h(0- f PQ{^,y)dC{y) d^ dx = [ ip{^) h{^) - ( Pp^{Yy) df{y) 

Jn Un L Jan J J in L Jan 

and so (34) holds a.e. and in fact everywhere thanks to Lemma 19 below. □ 

Lemma 19. Take a > 0 such that 2s -I- a ^ N and f G G[^^(fl). If u G L^{H,6{x)dx) solves 

{-^\nT'^ = f inVfn), 

then u G C'j^j’J^“(f2), the above equation holds pointwise, and given any compact sets K CC K' CC O, 
there exists a constant C = C{s,N,a,K,K',Tl) such that 

||M||c2«+a(K) < G (||/||c“(iC') + l|M||Li(n,5(x)tia;)) • 
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Similarly, if f £ and a € (0,2s), 

||M||c“(if) < C {\\f\\L’^(K') + \\u\\L^n,s{x)dx)) ■ 

In particular, if h is s-harmonic, then h £ C°“(f2) and the equality (—A|Q)*/i(a;) = 0 holds at every 
point X £ fl. 

Proof. We only prove the former inequality, the proof of the latter follows mutatis mutandis. Given 
X £ let 

v(x) = [ G\(f\x,y)u(y) dy. 

Jn 

Observe that v is well-defined and 


II L^{Q,,5{x)dx) < C'(^)-^1 s)ll'«||Li(n.5(x)da:) (36) 

Indeed, letting ipi > 0 denote an eigenvector associated to the principal eigenvalue of the Laplace 
operator, it follows from the Fubini’s theorem and Lemma 9 that 

f Pi(x) f G}(f'‘(x,y)\u{y)\dydx = \l~^ ( \u(y)\tpi(y) dy. 

Jci JQ J Q, 

In addition, —An = / in T>'{H), since for ip £ C)P(fl), 


/ n(-A)</j= / u(-AI^)"(p, 
In Jn 


thanks to the Fubini’s theorem, equation (24), Lemma 9 and Definition 15. Observe now that if 
C'^(Vl), then 


(-Ay^-V = 


^+oo 


r(i-s) 




— p. 




dt. 


ip £ 


The above identity is straightforward if ip is an eigenfunction and remains true for (p £ C^(Il) by density, 
using the fast decay of spectral coefficients, see (22). So, 


/ uip dx = / n(—A| q)^ ‘‘ipdx = 
Jq Jq, 


r(l — s) Jn Jo 


vP 


-1 P-e-' 


up 


t 


and 


Choose f £ C(: 


^+oo 


U = 


±s-l 


r(i-s) 


r(i-s) JnJo 

V — 


dt dx = 


ipP 


V — e 


dt dx 


t 


dt. 


such that / = / in K', ||/||c-(r«) < C'||/||c“(if') and let 
u{x) = CN,sf \x-y\~^^~^'‘^f(y) dy 


solve (—A)®u = / in K.^. It is well-known (see e.g. [19]) that ||M||c 2 s+<»(RiV) < C'||/||c£«(rn), for a constant 
C depending only on s, a, N and the measure of the support of /. It remains to estimate u — u. Letting 


7 {x)=cn,i-s f \x-y\ ‘'^^u(y)dy, 


we have as previously that —Av = / and 

_ s 

u = 


n-\-oo 


j.S-1 


r(i-s) 


V — e 


dt, 


10 












where this time e *'^v{x) = /rw e « v{y) dy. Hence, 

dt 




'0 


t 


^+oo 


+ / r 

Jo 




c?t. 


Fix a compact set K" such that K CC K" CC K'. Since i; — t; is harmonic in K', 
Ik — ^llc 2 «+“+ 2 (if") < C\\v — v\\l^(^k') < C\\u\\L^(^Q^s{x)dx) 


By parabolic regularity, 


r-‘( 


(v — v) — e — v) 

t 


dt 


020 + “ (iC) 


— ^Ik II (Q,5(3;)rfa:) ■ 


In addition, the function w = — e~*^v solves the heat equation inside fl with initial condition 

r(;(0, •) =v — v. Since A(r;—U) = 0 in iF', it follows from parabolic regularity again that r(;(t, x)/t remains 
bounded in (7^®+“ (it') as t —0+, so that again 


+ 00 


j.s-1 


3 — P **^7 


dt 


< c|kl 


C^^+°‘{K) 


L^{^^5{x)dx) ■ 


□ 

Lemma 20. Take a > 0, a ^ N, and u € n L^{td,S{x)dx). Given any eompact set K CC 

K' CC H, there exists a constant C = C{s, N,a, K, K',n) such that 

ll(“^ln)'''“llc“(if) < C* (Ikllc 2 -+°(_fs'') + IklkhS^.'Sla^ida:)) ■ 


Proof. With a slight abuse of notation, we write 

= f G]ff‘'{x,y)u{y)dy. 

Jn 

By Lemma 19 we have 

ll(“^la)'* ^'“Ilc2+“(if) < A (^|kllc2«+<»(if') + IK-Ak)® ^M||Li(n,5(x)dx)) < 

< G (Ikllc2s+“(ic') 

Obviously it holds also 


klkhn,<5(a;)dx)) ■ 


By (24), 


(-A)o(-Ak)^-^«||c.(x) < ||(-Ay*-^zi|k.+.(K). 
(-A) o{-A\^r-\ = {-A\^ru in V'{^\ 


which concludes the proof. 

Proposition 21. Let f € L^{id,S{x)dx) and u £ L\^JLl). The function 


□ 


u{x) 


Gtiix,y)f{y) dy 


belongs to LP(ri,S{x)dx) for any p C 


1 , 


A^ + 1 \ 

A + 1 - 2s 7 ■ 
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Proof. We start by applying the Jensen’s Inequality, 


G'a(x, 2 /)/(j/) dy 


< 


II {Q.S{x)dx) 


[ 

Gaix,y) 

In 

Hy) 


S{y)f{y) dy, 


so that 


|u(a:)|^(5(a;) dx < 


lLi(n.5(x)da:) / 
yenJn 


G!i{x,y) 

s{y) 


p 

d{x) dx 


and by (30) we have to estimate 


sup 

yen 



1 

y^(N- 2 s)p 




Pick £r > 0. Clearly, 


sup 

{y : 5 (y)>£} 



1 

y^{N-2B)p 


S{x)P+^ 

[\x-y\2 + S{x)6{yr 


dx < Cg. 


Thanks to Lemma 38, we may now reduce to the case where the boundary is flat, i.e. when in a 
neighbourhood A of a given point y € such that S{y) < e, there holds A n dO. C = 0} and 
A n O C {yjv > 0}. Without loss of generality, we assume that y = (0,?/Ar) and x = {x', xn) G 
B X (0, 1) C X R. We are left with proving that 


Ib Io [|a;'p + kiv - 

is a bounded quantity. Make the change of variables xn = 
= yNP- Then, the above integral becomes 


j+i 


[\x'\2 + \XN - ywp + XNyNf 
y]\[t and pass to polar coordinates in x', with 


-(Af+l-2s)p+Ar+l 

dN 


r^/VN r 

/ 

JO Jo 


1/yw 


„N -2 


dt 




[p2 + _ ^2-^{N-2b)p/2 [p2 _ 1|2 + 


P ■ 


(37) 


Now, we split the integral in the t variable into Jq'^^ + ■ Note that the exponent —{N + 

1 — 2s)p + TV + 1 is positive for p < (TV + 1)/(TV + 1 — 2s). We drop multiplicative constants in the 
computations that follow. The first integral is bounded above by a constant multiple of 


A/vn 


r-i/yw 


dp 


r.l/2 


„N -2 


dt 


fP+i 


0 [p2 + l]W-2«)p/2 [p2 + l+tf 


< 


A/vn 

Jo ' 


r.N-2 


dp 


[P^ + 1] 


(Ar+2-2s)p/2 


which remains bounded as y^v \- 0 since 

TV- 1 

p > 1 > 


TV + 2 - 2s 
The second integral is of the order of 


implies (TV + 2 — 2s)p — TV + 2 > 1. 


A/vn 


r3/2 


„N -2 


< 


dp 

AIvn 

3 

AIvn 

\ 

3 

r‘C30 

dp 


dt 


1/2 [p2 + |t_l|2](^-2*)p/2 [p2 + lf 


dp 


r.l/2 


„Af-2 


dt 


0 [p2+i2]W-2*)p/2 [p2 + lf 


P 


,Ar-l-(Ar-2s)p Al^’^.p) 


[p2 + 1]^ 

N —1 — {N—2s)p 


dt 


1 


[l+t2] 


(N-2s)p/2 


[p2 + 1]^ 
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which is finite since p < {N + 1)/{N + 1 — 2s) < N/{N — 2s) implies — 1 — (N — 2s)p > —1 and 
p > 1 > TV/(TV + 2 — 2s) implies 2p — TV + 1 + (TV — 2s)p > 1. 

We are left with the third integral which is controlled by 


r^tvN pi/yN tP+i 

lo [p2+i2]W+2-2.)p/2 - 


< 


N — {N-\-l — 2s)p 


r^/vN 

< dp p 

Jo 

r'^/vN 

< dp pN-{N+l-2s)p 

Jo 


I 


i/iVNP) 


dt 


fp+1 


3/(2p) 


[l+t2] 


(Af+2-2s)p/2 


The exponent TV — (TV + 1 — 2s)p > —1 since p < (TV + 1)/(TV + 1 — 2s), so this third integral is bounded 
above by a constant multiple of y“i“^+(^+i“ 2 s)p simplifies with the factor in front of (37). □ 


4 Boundary behaviour 

We first provide the boundary behaviour of the reference function hi. Afterwards, in Proposition 23 
below, we will deal with the weighted trace left on the boundary by harmonic functions induced by 
continuous boundary data. 

Lemma 22. Let hi be given by (8). There exists a constant C = C(TV, n,s) > 0 such that 

1 ^-(2-2s) pjg-{2-2,)^ (33) 

o 

Proof. Restrict without loss of generality to the case where x lies in a neighbourhood of dLl. Take 
X* € dft such that |a: — a:*| = S{x), which exists by compactness of dil. Take T C dLl a neighbourhood 
of X* in the topology of dLt. By Lemma 38 in the Appendix, we can think of T C {xn = 0}, x* = 0 and 
X = (0, 5(x)) G R-^~^ X K without loss of generality, in such a way that it is possible to compute 

f _ M _ dP 

'sJ|z'|2 + 5(a;)2]Af/2+i-. ■ 

Recalling (25), we have reduced the estimate to 


, _ ^\N+2-2s 



da{z) X 


X [ __ dz' X 

Jb^ [I^T + (5(x)2]^/2+l-- 


r 5{x)t^-^ 

Jo [t2 + 5(a;)2]iV/2+l-s 


r/5{x) 


5{x)^t 


N^.N-2 


[(5(x)2t2 + 5(^)2]iV/2+l- 


dt 


S{x) 


2s-2 


r/6{x) 


±N-2 


dt 


[t2 + lJJV/2+1- 


and this concludes the proof, since 


f'T / 5{x) 


j.N-2 


dt 


[t2 + 1]W/ 


2+1-s 


1 . 


□ 


In the following we will use the notation 

Pol? := [ P^{;0)gie)da{9) 
Jon 

where a denotes the Hausdorff measure on dLt, whenever g G L^(fl). 
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Proposition 23. Let ( € C(dfi). Then, for any z G dVl, 


^C(a^) xgn 


hi{x) 


> C(z) uniformly on dVL. 


(39) 


Proof. Let us write 
PoC(a:) 




-C(^) 




< 


1 


hi{x) 


Pi 


'an 


( PMx,e)m 

I an 

,ix,e)\aO)-az)\da{9)<CS{xf-^^ [ 

Ja 


an F 




_ Q\N+2-2s 


< C6{x) [ 

Jan 


\m-az)\ .... 

\x-0\^ 


It suffices now to repeat the computations in [1, Lemma 3.1.5] to show that the obtained quantity 
converges to 0 as x —>■ z. □ 

With an approximation argument started from the last Proposition, we can deal with a ^ G L^{dLl) 
datum. 

Theorem 24. For any f € L^(dfl) and any </> G it holds 

^ J{s{x)<t} h\{x) ito Jan 

Proof. For a general ( G L^{dLl), consider a sequence {CfcjfcgN C{dLl) such that 

/ \Ck{y)-Ciy)\d(T{y)-^0. 

Jan 

For any fixed /c G N, we have 

' ^ , , 4>{x) dx — / (j){x) C,{x) da{x) 

T^iix) Jan 

FUi^)-Ff,Ck{x) 


(40) 


t J{S{x)<t} 


< 


{Six)<t} hi{x) 

1 r ^fiUx) 


4>{x) dx 


(j){x) dx — 4>{x) Ck{x) da{x) 

Jan 


t J{ 5 {x)<t} di{x) 

+ / (j){x) Ck{x) da{x) - / (j){x)C{x) da{x) 

J J 

Call Afc := (fk — C- the term (41) equals 


(41) 

(42) 

(43) 


1 r P^Afc(x) 

t J{5{x)<t} hi{x) 


(j){x) da{x) = 


I an 


Call 


^(t.y) := 7 / 
t J(S 


1 r P^ix,y) 

t J{5{x)<t} hi{x) 


cj){x) dx Xk{y) da{y) 


/{S{x)<t} hi{x) 

Combining equations (25), (38) and the boundedness of (p, we can prove that 4) is uniformly bounded in 
t and y. Indeed, 


l^(t,y)l < 


ll<(>IU“(a) 


t J{5{x)<t} \^ - y\ 


s{x)^ ^ ll<(>IU”(a) 


N-\-2-2s — 


t J{s(x)<t} - y\ 


d{x) 

N 


dx 
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and reducing our attention to the flat case (see Lemma 38 in the Appendix for the complete justification) 
we estimate (the ’ superscript denotes an object living in 


XN 


0 Jb' [|x'|2 + X%] 


N/2 


dx' dxN < — ( I 

^ Jo Jb' , 


dC 


[I^P + 1] 


N/2 


dxN < 


Thus ^{t,y)Xk{y)da{y) is arbitrarily small in k in view of (40). 
The term (42) converges to 0 as t j, 0 because the convergence 


JrN -1 [|^|2 + i ]^^/ 2 - 


^hCkjx) 

hi{x) 




^ Ck{z)(l>iz) 

—¥Z 


is uniform in z S dfl in view of Proposition 23. 

Finally, the term (43) is arbitrarily small with k t +cxd, because of (40). This concludes the proof of 
the theorem, because 


lim 

t-i-O 


(l){x) dx — [ (j){x) C,{x) da{x) 
i(a^) Jan 

< ll^llz,“((o,to)xan) / lCfc(2/)-C(y)l c?cr(y) + ll<i^lU“(an) / \Ck{y) - C,{y)\d(j{y) 
J dci j on 


'{(5(x)<t} 


and letting k t +oo we deduce the thesis as a consequence of 40. 
Moreover we have also 


□ 


Theorem 25. For any y G A4{id), such that 


S d\y\ < oo, 


(44) 


and any (j) G C'*^(f2) it holds 


'{5(x)<i} 


hi{x) 


4>{,x) dx 


t^O 


(45) 


Proof. By using the Jordan decomposition oi y = yF — y~ into its positive and negative part, we can 
suppose without loss of generality that y > 0. Fix some s' G (0,s A 1/2). Exchanging the order of 
integration we claim that 


'{S{x)<t} 


Gh.{x,y)5(x) "dx^ 


Cf-^^'Siyf^' 5{y)>t 


CtS{y) 


Siy) < t 


(46) 


where C = C{N^ s) and does not depend on which yields 


^ Jo, \J{s{x)<.ty 


GUx,y)-^]dy{y) < C[ S{yf^'dy{y) 
hl[X) J J{S{y)>t} 


C / 5{y) dy{y). 

J {(5(a:)<i} 


The second addend converges to 0 as t j, 0 by (44). Since <^( 2 /)^® converges pointwisely to 0 in 
n as t j, 0 and < 5{y) in {S{y) > t}, then the first addend converges to 0 by dominated 

convergence. This suffices to deduce our thesis (45). 

Let us turn now to the proof of the claimed estimate (46) . For the first part we refer to [11, Proposition 
7] to say 


[ Gl,{x,y)5{xf-^^dx<t^-^^' [ Gl,{x,y)5[x)-^^+^^'dx<Ge-^^'5^^'. 

'J {S{x)<t} J{6(x)<.t\ 




We focus our attention on the case where dfl is locally flat, i.e. we suppose that in a neighbourhood 
A of y it holds A fl did C {xn =0} (see Lemma 38 in the Appendix to reduce the general case to this 
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one). So, since 6{x) = xm and retrieving estimate (30) on the Green function, we are dealing with (the 
' superscript denotes objects that live in R^“^) 


lo Ib'(v') W - y' 


3-2s 
VN Xj^ 


dx' 


{y') W - y'V + {XN - VnY + Xn yN [|a;' - y'P + {XN - 


dXN- 


From now on we drop multiplicative constants depending only on N and s. Suppose without loss of 
generality y' = 0. Set xn = yNV and switch to polar coordinates in the x' variable: 


rt/vN rl 




r,N-2 


dr 


Jo Jo r^ + y%{r]-l)^ + 'ny% [^,2 _ i) 2 ](^- 2«)/2 

and then set r = y^p to get 


drj 


rt/vN 

y% / 

Jo 


r'i-ivN 


r.N-2 


dp 


0 [p2 + (,y - 1)2](^-2-)/2 p2 + _ 1)2 + ^ 


dp < 


dp 


i-tlvN c'^/yN ^ 

< 7/2 / 773-2* / _ P _ 

-^^70 ^ Jo [p2 + _ 1)2](3-2.)/2 ^2 + (^_1)2+^ 


dp. 


Consider now s G (1/2,1). The integral in the p variable is less than 

r^/VN p 

SO that, integrating in the p variable, 

rt/vN 


— 3-|-2s 


yh 


r^/VN 

/ p^-^^\p-l\-^+^^ dp<tyN 

Jo 


and we prove (46) in the case s G (1/2,1). Now we study the case s G (0,1/2]. Split the integration in 
the p variable into and the latter can be treated in the same way as above. For the other one 

we exploit the inequality + (?7 — 1)^ + ?7 > | to deduce^: 


2/w / n 

Jo 


3-2s 


r-l/yiv 


dp 


lo [p 2 + (77-l)2+7y](3-2*)/2 p 2 + _ 1)2 

i-'i./yN 

< Vn I P' ' 


dr] < 


3-2s 


lo [p2 + (,y-l)2] 


dv<y%J^ 


2 p3-2s 


\v-M 


l-2s 


dp < Pn- 


Note now that, in our set of assumptions, y^ = 5{y) < t. So 7 /^ < ty^ and we get to the desired 
conclusion (46) also in the case s G (0,1/2]. □ 


5 The Dirichlet problem 

Recall the definition of test functions (7). 

Lemma 26. T(n) C (70(12) (7 (7°“(r2). Moreover, for any ip G T(12) and z G did, 

= J^Pn(.y,z){-A\^y^{y) dy. (47) 

^In the computation that follows, in the particular case s = ^ the term |1 — 77 ^®“^ must be replaced by — In |1 — 77 I, but 
this is harmless. 
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Proof. Take if € T(0) and let / = (—A\^Yif. Since / G CfP{n), the spectral coefficients of / have fast 
decay (see (22)) and so the same holds true for if. It follows that if G and T(n) Q CQ{fl). By 

Lemma 9, for all x G 

if{x) = [ Gliix.y) f{y)dy 
Jn 

Using Lemma 12, (28) and the dominated convergence theorem, (47) follows. 

Since (—Alj.^)® is self-adjoint in H{2s), we know that the equality (—= / holds in T>'(U) and 
the interior regularity follows from Lemma 19. □ 

Lemma 27 (Maximum principle for classical solutions). Let u G fl L^(U, 5{x) dx) such that 

(—A|q)*u>0 in LI, liminfM(a;) > 0. 

x—^dQ 

Then u > 0 in LI. In particular this holds when u G T{Ll). 

Proof. Suppose x* G LI such that u{x*) = < 0- Then 

= [ [u{x*)-u{y)]J{x,y) dy + k{x*)u{x*) <0, 

Jn 


a contradiction. □ 

Lemma 28 (Maximum principle for weak solutions). Let y G Ai{Ll), ( G A4(dLl) be two Radon measures 
satisfying (9) with y, > 0 and C > 0. Consider u G Lj^Y^) ® weak solution to the Dirichlet problem (10). 
Then u > 0 a.e. in LI. 

Proof. Take / G CY’(Ll), / > 0 and if = {—A\^)~^f G T{Lt). By Lemma 27, if > 0 in Lt and by Lemma 
26 —^ > 0 on dLl. Thus, by (11), J^uf > 0. Since this is true for every / G C^{Ll), the result 
follows. □ 

5.1 Proof of Theorem 5 

Uniqueness is a direct consequence of the comparison principle. Lemma 28. Let us prove that formula 
(12) defines the desired weak solution. Observe that if u is given by (12), then u G L^{Ll, 5{x)dx). Indeed, 


/ Tii.x) / 
J Q J Vt 


GQ{x,y)dy{y) 


dx < 


Cl Jo, 


Gn{x,y)(pi{x)dx d\y\{y) 


= ^ [ Ti{y) d\y\{y) <G\\Sy\\M{<n) (48) 
"^1 Jn 


This, along with Lemma 18, proves that u G L^(Ll, d(x)dx) and (13). Now, pick if G T{Ll) and compute, 
via the Fubini’s Theorem, Lemma 9 and Lemma 26, 


«(a;)(-^ln)"V'(a;) dx = 

= f if GUx,y)dyiy)] i-A\^Yifix) dx+ f (f P^{x, z) dC{z)] {-A\^yif{x) dx = 

= f if Ghix,y){-A\^Yif{x)dx'\dy{y)+ f ff P^{x, z){-A\^Yif{x) dx^ dC{z) = 

Jn \Jn / Jon \Jn / 

= [ Y{y)dy{y)- [ ^{z)dC{z). 
Jn Jon ov 

□ 
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5.2 Proof of Lemma 4 

Proof of 1. Consider a sequence {?7fc}feeN C of bump functions such that 0 < r]i < ... < rjk < 

Vk+i < • ■ • < 1 and r]k{x) f Xn{x) as k f oo. Consider G aird define fk '■= € 

Yk-.= i-A\^rfk€Tin). 

Let us first note that the integral 

Jn 

makes sense in view of (31) aird (13). The sequeirce {fk}k^fi trivially coirverges a.e. to (—A|q)*' 0, while 


\Yk{x) - Y{x)\ < [ G^(x,2/)|(-A|s.j)X2/)|-(l-77fc(2/))dj/ 

Jn 

converges to 0 for any a; G 17 by dominated convergence. Since it is a weak solution, it holds 

[ ufk = [ 'tpkdfx- f dC. 

Jn 7a Jan 

By dominated convergence JqU fk Jq u( “ A|j.j)®i/) and ipk ^ /n Lrdeed for any fc G N, 


|/fc| < |(-A|fj)>| and \Yk\< 


Af 


(-A|n)V 




Pi. 


L°°{Q) 


where the latter inequality follows from the maximum principle or the representation formula ijjkix) = 
InGQ{x,y)fk{y)dy. Finally, the convergence 


f dfik 

Jan 


dC 


fc't'OO 


0 


holds by dominated convergence, since 


f dipk 

Ian dv 


dc= / mfk 

Jn 


by the Fubini’s Theorem, and PqC G L^{n,6{x)dx) while \fk\ < |(—Alj.^)®^! < C6 by Lemma (16), so 
that 

/ mfk^ [ mf = 0 

Jn ^=1°“ Jn 

because is s-harmonic and / G 

The proof of the boundary trace can be found in Theorems 24 and 25, by recalling the representation 
formula provided by Theorem 5 for the solution to (10). □ 

Proof of 2. Recall that u is represented by 


u{x)= f Gf^{x,y) y{y) dy + f P^{x,y) C,{y) da{y). 

Jn Jan 

By Point 1. and Lemma 19, u G Moreover, u G L^{Pt,5{x)dx) thanks to (13). So, we can 

pointwisely compute ( —A|Q)®it by using (3) and (54): this entails by the self-adjointness of the operator 
in (54) that 

/ (-A|n)'’wV’ = f u(-A|s.j)>= [ yif. for any i/-G G“(17) 

Jn Jn Jn 

and we must conclude that (— Aj^)*!! = y, a.e. By continuity the equality holds everywhere. 

We tmir now to the boundary trace. The coirtribution giveir by Gq/t is irrelevant, because it is a 
bounded fmrctioir as it follows froirr 


|G^m(^)I < G||mIU==(0) / -— % 

Jn \x — y\ 
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where we have used (30). Therefore, by Propositions 23, there also holds for all xg G dfl, 


/2l[X) 


lim 

x—>xo,xG^ 


+ ^hC(x) 

hi(x) 


= C(xo)- 


□ 

Proof of 3. By Lemma 20 /i S In addition, we have assumed that ^ £ (7(912). Consider 

v{x) = [ Gt^{x,y)^i{y) dy+ [ P^{x, z) C(,z) da{z) 

Jn Jon 

the weak solution associated to data y and C,. By the previous point of the Lemma, n is a classical 
solution to the equation, so that in a pointwise sense it holds 

u — V 

(—A|q)'*(m — n) = 0 in 12, —— = 0 on 912. 

hi 

By applying Lemma 27 we conclude that |m — n| < ehi for any e > 0 and thus u — v = 0. □ 


6 The nonlinear problem 

Lemma 29 (Kato’s Inequality). For f £ L^{id,S{x)dx) let w £ L^{n,6{x)dx) weakly solve 

ji-^\nYw = f inn 

I ^ = 0 on 912. 

I hi 

For any convex $ : R —^ K, $ £ (7^(]R) such that d>(0) = 0 and ^(ri;) £ L\^^{n), it holds 
Moreover, the same holds for $(i) =2+ = 2 A 0. 

Proof. Let us first assume that / £ (7;“^(12). In this case, by Lemma 19, w £ (7fj),’^“(12) and the equality 
(—A|j.j)®w = / holds in a pointwise sense. Then 


(— A|j.j)®$ o r(;(a:) = / [^{w{x)) — ^{w{y))] J{x,y) dy + k{x) ^{w{x)) 

Jn 

= ^'{w{x)) / [w{x) — w{y)\J{x,y) dy + k{x)^(w{x)) 

Jn 

- [ [w{x) -w{y)Y J{x,y) [ {w{x) + t[w{y) - w{x)]){l - t) dt dy 

Jn Jo 


< $'(ii;(a;)) (-A|s. 2 )"n;(x) 


where we have used that 3>" > 0 in R and that 3>'(2) < 2$(2), which follows from $(0) = 0. 

We deal now with / £ L°°(12). Pick {fj}j^^ C (7“(12) converging to / in L^(12, 6{x),dx) and bounded 
in L°°{n). The corresponding {wj = converges to w in L^{n,S{x)dx), is bounded in L°°(12) 

and without loss of generality we assume that fj —>■ / and Wj —>■ w a.e. in 12. We know that for any 

£ r(12), V' > 0 

JQ J Q 

By the continuity of $ and $' we have ^(uij) —2 $(ui), ^'(wj) —>■ ^'{w) a.e. in 12 and that {‘f’(n'j)}jg|!j, 
are bounded in L°°(12). Since is converging to / in L^{n,S{x)dx), then 


/ (-A|n )>—2 [ ^{w){-A\^Yfj and f 

Jn Jn Jn 


by dominated convergence. 


f<^'{w)tp 
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For a general / G L^{Q,S{x)dx) define fj^k ■= (/ A j) V {—k), j,k £ N. Also, split the expression of 
d) = <I)i — <1)2 into the difference of two increasing function: this can be done in the following way. The 
function <&' is continuous and increasing in R, so that it can either have constant sign or there exists 
to G R such that <i)'(to) = 0. If it has constant sign than <!> can be increasing or decreasing and we can 
choose respectively <I)i = <!>, <1)2 = 0 or <I)i = 0, <1)2 = —‘h- Otherwise we can take 


$i(t) 


$(t) t> to 
d>(to) t <to 


and < 1*2 (t) 


0 t > to 

d>(to) - «>(t) t<to' 


We already know that for any ip G T(0), ip > 0 


Jn JQ, 

On the right-hand side we can use twice the monotone convergence, letting j t oo first and then fc f oo. On 
the left hand side, by writing <1) = <l)i — < 1)2 again we can exploit several times the monotone convergence 
by splitting 


J J ^ J Q, 

- [ <i>2Kfc)[(-A|n)>]++ / <l)2Kfc)[(-A|„)>]- 
J Cl JQ 

to deduce the thesis. 

Finally, note that ^(t) = t~^ can be monotonically approximated by 





1 


which is convex, and <l)j(0) = 0. So 


/ (-A|f2)*V'< / f^^'jMip. 

J Q JQ 

Since ^j{t) 1 1^ and 2<1)'(t) t 1 + sgn(t) = 2x(o, ,+oo)(i)) we prove the last statement of the Lemma. □ 

Theorem 30. Let f{x,t) : O x R —> R 6e a Caratheodory function. Assume that there exists a 
subsolution and a supersolution u,u G 6(x)dx) H to 


(— A|j.;)®tt = /(x, u) in n 
u 

— =0 on dLl 

hi 


(49) 


Assume in addition that f{-,v) G L^{Ll,5(x)dx) for every v G L^{fl,6{x)dx) such that u < v <u a.e. 
Then, there exist weak solutions ui,U 2 G L^{n, 5(x)dx) in \u,u] such that any solution in the interval 
[u, m] satisfies 


U<U 1 <U<U 2 <U a.e. 

Moreover, if the nonlinearity f is decreasing in the second variable, then the solution is unique. 
Proof. According to Montenegro and Ponce [17], the mapping v F(-,v), where 

F(x, t) := f{x, [t A ?l(x)] V u{x)), x G O, t G R, 
acts continuously from L^{Vl,5{x)dx) into itself. In addition, the operator 


S{x)dx) 

JC{v){x) = ( Gli{x,y)F{y,v{y))dy 
JQ. 


1C: L\n,S{x)dx) 
v(x) 


20 



is compact. Indeed, take a bounded sequence {'I’nljjgpj in L^{i^,S{x)dx). On a compact set K CC O, 
u,u are essentially bounded and so must be the sequence By Theorem 5 and Lemma 

19, bounded in 0 LP{n,6(x)dx), p G [1, (iV + l)/(iV + 1 — 2s)). In particular, a 

subsequence converges locally uniformly to some v. By Holder’s inequality, we also have 

Ikrifc - v\\L^n\K,Six)dx) < WVrik “ II LP(n\if ,5(a:)da;) II ln\iC II Lp'( n\if ,5(a:)da;)' 

Hence, 

Ikrifc - v\\L'-ISl,5{^)dx) < Iknfc - v\\(^K,5(x)dx)\\S\\L^{n) + C'll ln\iC IIlp'( n\if,5(a,)dx' 

Letting k —>• +oo and then K ^ fl, we deduce that /C is compact and by the Schauder’s Fixed Point 
Theorem, 1C has a fixed point u G 5(x)dx). We then may prove that u < u < u hy means of the 

Kato’s Inequality (Lemma 29) as it is done in [17], which yields that u is a solution of (49). 

The proof of the existence of the minimal and a maximal solution ui,U 2 G 6(x)dx) can be 

performed in an analogous way as in [17], as the only needed tool is the Kato’s Inequality. 

As for the uniqueness, suppose / is decreasing in the second variable and consider two solutions 
u,v G L^{d,5{x)dx) to (49). By the Kato’s Inequality Lemma 29, we have 

(-^ln)*('“ - ^ X{u>v}[fix,u) - fix,v)] <0 in H 

which implies {u — v)~^ < 0 by the Maximum Principle Lemma 28. Reversing the roles of u and v, we 
get also {v — m)+ < 0, thus u = v in Cl. □ 


6.1 Proof of Theorem 6 

Problem (17) is equivalent to 

J (PnC -v) in H 

I ^ = 0 on 

I hi 


(50) 


that possesses u = PqC as a supersolution and u = 0 as a subsolution. Indeed, by equation (38) we have 

0 < P?^c < IICIUp»(Q)/ii < C|[C|[i=.(a)^-("-'*^ 

Thus any v G L^{Cl,6{x)dx) such that 0 < u < ^nC satisfies 


g{x,v) < h{v) < h{cS G L^{Cl,6{x)dx). 

So, all hypotheses of Theorem 30 are satisfied and the result follows. □ 


7 Large solutions 

Consider the sequence built by solving 

(-A[^)% = -u^ 

Ml 

h=^ 
hi 

Theorem 6 guarantees the existence of such a sequence if G L^{Cl,5{x)dx), i.e. p < 1/(1 — s). 

We claim that the sequence is increasing in Cl: indeed the solution to problem (51) is a subsolution for 
the same problem with boundary datum j + 1. In view of this, the sequence {ujljgN admits a pointwise 
limit, possibly infinite. 


in Cl 
on dCl. 




(51) 
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7.1 Construction of a supersolution 

Lemma 31. There exist So,C > 0 such that 


> -CS-°^P, forS<Soanda = 


2s 

p-1' 


Proof. We use the expression in equation (3). Obviously, 

= / [5(x)~°‘- S{y)~°‘]J{x,y) dy+ S{x)~°'k(x) > f [6(x)~°‘ - S{y)~°‘]J{x,y) dy. 

Jfi Jfi 

For any fixed x € ft close to dfl, split the domain ft into three parts: 

Oi = jy e O : S{y) > ^ 5(x)| , 

= jy e O : i S{x) < 6{y) < ^ S{x) 

^^3 = jy € O : S{y) < i (5(x)| . 

For y € Oi, since S{y) > S(x), it holds S(x)~°' — S(y)~°‘ > 0 and we can drop the integral on Oi. Also, 
since it holds by equation (33) 

w N C 

J[x,y) < 


- y\ 


N+2s ’ 


\-Oip 


the integration on 1^2 can be performed as in [2, Second step in Proposition 6 ] providing 
PV [ [^(a:)-“ - 5(y)-“] Jix,y) dy > -C= -CS{xy 

Jn2 

To integrate on O 3 we exploit once again (33) under the form 

S{x)S{y) 


J{x,y) <C ■ 


\x - y\ 


iV+2+2s 


to deduce 


f [(5(a;) °‘-S{y) °‘] J{x,y) dy > -C 6{x) f 
J fio J r 2 j 


5{y? 


/03 \x - y\ 

Again, a direct computation as in [2, Third step in Prooposition 6 ] yields 


iV+2+2s • 


[ [5(x)-“ - (5(y)-“] J{x,y) dy > -C(5(:r)-“-2« = -CS{x)-^p. 
J Qq 


□ 


Lemma 32. If a function v € L^{ft, S(x)dx) satisfies 

{-^\qYx e (-A|s^)®n(a:) > -Cv(x)p, when S{x) < Sq, (52) 

for some C,So > 0, then there exists u G L^{ft,S{x)dx) such that 

(—AIqYu^x) > —u{x)P, throughout ft. (53) 

Proof. Let A := C^Pp~'^'> V 1 and flo = {a: S 11 : 6{x) < Jq}) then 

(-Ay'* (An) > -(Ay, inOo. 

Let also y := A||(— A|Q)®n||ioo(QyQg) and define u = yGfil + An. On u we have 

(-A|n)®M = /r +A(-A|fj)'*n > A|(-A|fj)'*n| + A(-A|s.j)'*n > -iT throughout ft. 

□ 


22 






Corollary 33. There exists a function u G L^(fl,S(x)dx) such that the inequality 

(—A|q)®m > —id’, in fl, 
holds in a pointwise sense. Moreover, u x 

Proof. Apply Lemma 32 with v = ( 5 “ 2 s/(p-i). corresponding u will be of the form 

□ 


7.2 Existence 

Lemma 34. For any j G N, the solution Uj to problem (51) satisfies the upper bound 

Uj < u, in fl, 


where u is provided by Corollary 33. 

Proof. Write Uj = jhi — Vj where 

J i = Uhi-Vjf in 

I ^ = 0 on dfl. 

(. hi 

and 0 < nj < jhi. Since (jTii — vjY G we dednce that Vj G for any a G (0,2s). By 

bootstrapping Vj G C°°{il) and, by Lemma 19, also Uj G C^(fl). This says that uj is a classical solution 
to problem (51). Now, we have that, by the boundary behaviour of u stated in Corollary 33, Uj < u 
close enough to 917 (depending on the value of j) and 


(-A|o)" (m-Mj) > u^-uP, in 17. 


Since —11^ G (7(17) and — vF = —oo, there exists a;o G 17 such that Uj{xo)P — u(xo)p = 

m =: max 3 ;gQ (uj{x)^ — u(x)^). If m > 0 then (—A|q)^(m — Uj)(xo) > m > 0: this is a contradiction, as 
Definition 3 implies. Thus m < 0 and Uj < u throughout 17. □ 


Theorem 35. For any p G 



there exists a function u G L^(17, 5(x)dx) solving 


■-A\^ru=-up 

= +00 


in 17 
on 917 


both in a distributional and pointwise sense. 

Proof. Consider the sequence provided by problem 51: it is increasing and locally bounded 

by Lemma 34, so it has a pointwise limit u < u, where u is the function provided by Corollary 33. 
Since p > 1 + s and u < then u G L^(fl,S(x)dx). Pick now ip G (7)?°(17), and recall that 

(5“^( — Alj.j)'’'^ G L°°{U): we have, by dominated convergence, 

/ ( u[-A\^Y'ii), f uPf)f uPy 

Jn Jtoo Jq Jq jtoo Jq 

so we deduce 

/ w(-A|o)> = - / uPy. 

Jn Jo. 

Note now that for any compact K CC 17, applying Lemma 19 we get for any a G (0, 2s) 

llc“(X) ^ A + Il'Wj ||Li(n,(5(s)dx)^ ^ A + ||w||Li(n,5(x)dx)^ 
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which means that is equibounded and equicontinuous in C{K). By the Ascoli-Arzela Theorem, 

its pointwise limit u will be in C{K) too. Now, since 

= -uP mV'{n), 


by bootstrapping the interior regularity in Lemma 19, we deduce u € So, its spectral fractional 

Laplacian is pointwise well-defined and the equation is satisfied in a pointwise sense. Also, 


. r “( 2 ^) 

limmi ■;——^ 

hi{x) 


K-s-on hi[x) 


= j 


and we obtain the desired boundary datum. 


□ 


A Appendix 

A.l Another representation for the spectral fractional Laplacian 

Lemma 36. For any u G F[ (2s) and almost every a; £ 12, there holds 


= PV [ [u{x) - u{y)]J{x,y) dy + k{x)u{x), 
JQ 


where J{x,y) and k{x) are given by (4). 

Proof. Assume that u = ipj is an eigenfunction of the Dirichlet Laplacian associated to the eigenvalue 
Xj. Then, (— = A®u, = J^pnft, ■,y) u(y) dy = e~*^^u and ior all x G it 


EiL^^.A\^Yu{x) = 



u{x)- / Pn{t,x,y)u{y) dy 
Jn 


dt 

fTTs 


= lim / [ Pn{t,x,y)[u{x) - u{y)] dy+ [ u{x)[l- [ pn{t,x,y) dy] 

Jo \ Jq ) 

= PV / [u{x) — u{y)]J{x,y) dy + k{x)u{x) 

Ja 


dt 

fT+2 


(55) 


By linearity, equality holds on the linear span of the eigenvectors. Now, if it £ iJ(2s), a sequence 
of functions belonging to that span converges to u in H{2s). In particular, (—A|j.j)®u„ converges to 
{—A\^Yu in L^(f2). Note also that for v £ L^(12), 


s 




dt 

p+I 


• v{x) dx 


+ 00 

xiukVk 

fc=i 


< lkllff(2s)lkllL2(a) 


so that we may also pass to the limit in LJ{Vl) when computing (54) along the sequence By 

the Fubini’s theorem, for almost every a; £ 12, all subsequent integrals are convergent and the identities 
remain valid. □ 


A.2 The reduction to the flat case 

In this paragraph we are going to justify the computation of the asymptotic behaviour of integrals of the 
type 

[ -y\) dy as 5(x)|o, 

J Ann 
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where A is a fixed neighbourhood of x with A fl dfl ^ 0, by just looking at 


[ dt f dy'F(s{x),t,'/WW+W-5{^Y 

The first thing to be proved is that 

\x - j/p X Ixo - 2/oP + |<^(a::) - S{y)\'^, 
where xq, yo s-re respectively the projections of x, y on dfl. 

Lemma 37. There exists e = e(fl) > 0 such that for any x € fl, x = xq + S(x)vS(xo), xq € dfl, with 
S(x) < e and any y G fl with S{y) < e and \yQ — xqI < e 

^ (ko - 2/oP + |^(a:) - S{y)\^) < \x - y\^ < ^ (|a;o - yoP + \S{x) - 6{y)\^) . 

Proof. Call Oe = {a: G n : S{x) < e}. Write a; = Sq + (5(a:)V(5(a;o), y = yo + ^(y)Vi5(yo)- Then 

\x - yp = \xo - yo? + |<5(a;) - S{y)? + S{yf\v6(xo) - vS{yo)? + 2[(5(a;) - ^(y)](a;o - yo, ^S{xq)) + 

+ 26{y){xo - yo, ^S{xo) - V(5(yo)) + 25(y)[(5(a;) - 5(y)](V(5(a;o), vS{xo) - V(5(yo))- (56) 

Since, for e > 0 small, S G C^’^(fle) and 

lvS(x)-vS(y)?<Hdl?a^,.^^Jx-y? 

(xo - yo, 'fS{xo)) = O{\xo - yo?) 

|(a:o - yo, V(5(a;o) - V(5(yo))| < ||<5||ci.i(n,)|a:o - Vo? 

|(5(a:) -(5(y)| = ||(5||ci.i(n^)|a: - y| 

|(Vi5(a:o), Vi5(a;o) - V(5(yo))| < ||^||ci.i(n,)|a;o - yo| 

The error term we obtained in (56) can be reabsorbed in the other ones by choosing £ > 0 small enough 
to have 

d{y)'^\^S{xo) - V(5(yo)P + 2|(5(a;) - S{y)\ • |(a:o - yo, V(5(a:o))| + 2S{y)\{xo - yo, Vi5(xo) - V(5(yo))| + 

+ 2(5(y)|(5(a:) - 5(y)| • |(Vi5(a;o), Vi5(xo) - V(5(yo))| < ^ (la^o - Vo? + 115(2;) - 5{y)?) . 


□ 

Lemma 38. Let F : (0, +c»)^ (0, +oo) be a continuous function, decreasing in the third variable, 

and = {a; G n : d{x) < e}, with e = s{fl) > 0 provided by the previous lemma. Consider x = 
Xo + (5(a;)Vi5(xo), xo G dfl, and the neighbourhood A of the point x, defined by A = {y G fie '. y = 
yo + d[y)s/5{yo), |a;o — yo| < £}• Then there exist constants 0 < ci < C 2 , ci = ci(fl),C 2 = C 2 {fl) such that 

I I + \t- i5(a;)p^ < y F((5(a;),(5(y), |a; - y|) dy < 

< C 2 y dt j dy' F(^5{x),t,ci\/\y’\^ + 1^ - <5(a;)p^ . 
where the ’ superscript denotes objects that live in 

Proof. By writing y = yo + 5(y)Vi5(yo), yo G dfl and using the Fubini’s Theorem, we can split the 
integration into the variables yo and t = S{y): 

f F{5{x),S{y),\x - y\) dy < f ( [ F{S{x),t,\x - yo - tvS{yo)\) dt\da{yo). 

J A aBj(xo)nan \Jo / 
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Using the monotony of F and the above lemma, we get 

[ - y\) dy < [ ( [ dt\da{yo) 

JA JB^{xa)nan \Jo ^ ' J 

where c is a universal constant. Representing i?£(xo) n9U via a diffeomorphism 7 with a ball R' C 
centered at 0, we can transform the integration in the j/q variable into the integration onto R'. The volume 
element IR 7 I will be bounded above and below by 

0 < Cl < IR 7 I < C2, 

in view of the smoothness assumptions on dfl. □ 
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